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The fluctuations in the temperature and polarization of the cosmic microwave background 
are described by a hierarchy of Boltzmann equations. In its integral form, this Boltzmann 
hierarchy can be converted from the usual Fourier-space base into a position-space and causal 
description. We show that probability densities for random flights play a key role in this 
description. The integral system can be treated as a perturbative series in the number of steps 
of the random flights, and the properties of random flight probabilities impose constraints 
on the domains of dependence. We show that, as a result of these domains, a Fourier-Bessel 
decomposition can be employed in order to calculate the random flight probability densities. 
We also illustrate how the H-theorem applies to the cosmic microwave background: by using 
analytical formulae for the asymptotic limits of these probability densities, we show that, as 
the photon distribution approaches a state of equilibrium, both the temperature anisotropics 
and the net polarization must vanish. 
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I. INTRODUCTION 



The cosmic microwave background (CMB) has been a prolific source of theoretical [1| an d exper 



imental 



fl 



J] results over the last couple of decades. Because it relies on low-energy interactions well 
inside the linear regime, the CMB now occupies a central part in our understanding of Cosmology, 
providing crucial information about the initial state of the Universe and about many processes that 
distorted it since the surface of last scattering. 

A fundamental description of the temperature and polarization of the CMB starts with the 
basic interaction dynamics (free propagation and scatterings with matter sources), and is fully 
realized through the Einstein-Boltzmann hierarchy of equations. We shall here describe the system 
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in a backward fashion: starting with the integral version of Boltzmann's equations, we will uncover 
details about the physical system described by that set of equations. Our procedure will be directed 
by the underlying structure of those integral equations, in particular by the way in which the 
hierarchy couples the equations. By removing and then reintroducing the source terms we will be 
able to reveal some hidden details about the microscopic dynamics of the physical processes that 
lead to the properties of the CMB which are observed today. To be clearer, we shall demonstrate 
that the Boltzmann equations preserve, in a very fundamental way, the fact that the fluctuations in 
the temperature and polarization of the CMB were imprinted through Thomson scatterings of low 
energy photons and free electrons during the recombination. Consequently, the hierarchy can be 
rewritten in terms of a series expansion over the number of scatterings that photons suffered before 
decoupling from matter. In this series expansion, besides the probability for each interaction, the 
number of scatterings appears as a key ingredient, since it determines the probability density that 
mediates the way in which the source terms imprint their contribution into the final signal. 

The argument which leads to these conclusions will take us through a treatment of the CMB 
in position space and, in many senses, the present work is a further development of [§]. When 
performing the conversion from Fourier space to position space, it is unavoidable that a certain 
class of integrals over products of spherical Bessel functions appear in the description. The proper 
understanding and interpretation of these integrals constitute the key barrier that we should over- 
come in the description of the system in position space. As we shall see, these integrals codify the 
memory of the fundamental scattering dynamics. In fact, we shall demonstrate that these integrals 
constitute a generalization of random flight probability densities. 

Random flights are a classical problem in mathematical physics with many applications in 



physics and astronomy [3]. The problem was first ] 

in context of the study of bird migrations. Lord Rayleigh, soon after, applied the same ideas in 
acoustics. Further contributions on this subject are described in Ref. In very simple terms, 
a random flight (in a .D-dimensional space) is the trajectory performed by a body which moves 
at constant speed and changes its direction of motion into another random direction at Poisson- 
distributed time intervals. If the movement has a fixed origin, we may ask, based on the length of 
the intermediate displacements, as well as on the number of displacements, what is the probability 
for the moving body to reach a distance r from the origin. 

The precise form of the visibility function also plays a fundamental role in the description 
of the CMB in terms of random flights, since it models the probability of scatterings during the 
recombination - which changes as a function of time only. The spatial independence of the visibility 
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function and the elasticity of Thomson scattering (which are, of course, already codified in the 
usual Boltzmann hierarchy) are the fundamental ingredients which lead to a random-flight-based 
description of the CMB. 

This paper is organized as follows: after introducing the Boltzmann hierarchy in the integral 
form in Sec. [Ill we temporarily decouple the evolution of the temperature from the polarization, and 
show the consequences of this simplification for the expression of the temperature and polarization 
in Sec. IIIII We then show, in Sec. IIVI how the family of integrals over spherical Bessel functions 
that appears in that description is related to random flight probability densities functions. After 
clarifying the structure of this dynamical system under the approximation of uncoupled temperature 
evolution, we go back to the original form of the hierarchy in Sec. IVIH and show that the formalism 
introduced in the previous Sections can be extended to treat the general case. We present, in Sec. 
IVIII1 a computational tool based on Fourier-Bessel expansions which is suitable for computing the 
probability functions appearing in our description. Finally, in Sec. HXI we use asymptotic arguments 
to show that high order terms should not contribute to the observables, which is expected from 
Boltzmann's H-theorem. 



II. BOLTZMANN'S EQUATIONS 

The hierarchy of Boltzmann's equations describing CMB temperature fluctuations and polar- 
ization can be written in terms of a set of coupled integral equations - see, e.g., Q|. Let the 
temperature anisotropics observed at position x G and at (conformal) time r\ along the direction 6 
be given in terms of its momenta as: 

/ m 

Similarly, the polarization, in terms of the usual Stokes parameters Q, U and /, is decomposed as: 
Q + iU. _ f d 3 k 



G(x ,r ?0 ,6) = J ^0^^ 4n ^ l &i(KVo)YUk)Y lm (o) • (1) 



/d 6 k - 
— — e* kx ° 4vr i l «Kk, f]o) Yf m (k) 2 Y lm (6) , (2) 

^ ' Im 

where 2^lm are the spin- weighted spherical harmonics flo| . 
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The momenta of the CMB temperature and polarization are then given by the integral equations: 



@i(Kvo) 



>lo 



&sw(K rj)ji(k Ai] ) - V b (k, rj) jl(kArio) 



9 2 (k, rj) - VEa 2 (k,r}) 



+ + &)(k,r l )j l (kAr ]0 ) 



|if (kA m ) + \ji{kA m ) 



(3) 



and 



oci(k,rj ) 



3 (I 



(I -2)1 Jo 



*° dr] fi'^e-^) I [ @2 (k, n) - V6a 2 (k, rj) 



3l(kAr] 



(4) 



In Eqs. ©-([J]) a prime denotes a derivative with respect to conformal time rj, the optical depth 
to Thomson scattering is /j(r/), and we have defined the interval Atjq = rj — rj. It is sometimes 
convenient to define the ubiquitous source term P(k, rj) = ^ [02(k, r/) — \/6a2(k, rj)\ . The system 
above is closed once the perturbed Einstein equations are used to determine how the linear scalar 
cosmological perturbations 9sWi H> $ an d ^ evolve with time. However, both the precise nature of 
the perturbed Einstein equations, or of the initial conditions that are used to evolve those equations, 
are irrelevant for our results. 

Eqs. (J3J)— show that the primary sources of temperature fluctuations are the Sachs- Wolfe term, 
0SW, the baryon velocity, Vb, and the gravitational potentials $ and ^ (also known as the Bardeen 
potentials - we work in the conformal-Newtonian gauge). The primary source of the polarization 
of the CMB, on the other hand, is the quadrupole of the temperature fluctuations. The integral 
equations then couple all the momenta of temperature and polarization, mediated by the visibility 
function g(/j,) = /j / (r/)e~ /x ^^ . 

Henceforth we will take x D = 0, i. e., the observer is taken to be at the origin of the coordinate 
system employed for the description of the problem. 



III. UNCOUPLING THE TEMPERATURE EVOLUTION 

We can decouple Eqs. (EJ- ([4]) by neglecting the term a 2 in ([3]). This truncation represents the 
approximation whereby deviations from the equilibrium temperature are described by the Sachs- 
Wolfe (SW) and integrated Sachs- Wolfe (ISW) effects. Within this approximation, then Eq. ([3|) 
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becomes: 



6,(M ) = r° dr l {g(n)[Osw(Kv)ji(kA m )-V b O<:,r l )j' l (kA m )] 



+ 



(5) 



Neglecting polarization as a source term for the temperature anisotropies is in fact a very good 
approximation, and the reason for this underlies the argument presented in this paper. The visibility 
function g(rj) should be regarded as the probability per unit (conformal) time that photons will 
scatter with some free electron - and, in fact, g(r]) is defined in such a way that this probability 
is normalized, f °° drj g(rj) = J °°(i/ue _ ^ = 1. This means that each time a factor of the visibility 
function intermediates a source term, that source term is damped by a factor e, with < e < 1. Since 
the lowest-order polarization term has at least one factor of the visibility function, it contributes 
as a source term to the temperature with two factors of the visibility function. Hence, the SW and 
ISW effects dominate the intensity of the signal, and Eq. ([5]) accounts for the largest contribution 
to the temperature anisotropies. 

A. CMB temperature in position space 

Let's now take the lowest-order contribution to the temperature anisotropies, and express it in 
terms of position space. Expressing Eq. ([I]) as: 



e(7 ?0 ,6)=^^ ) (r ?0 )Y^(6), 



(6) 



Im 



the coefficients 0$(t] o ) are given by: 



x \v'(v) 



+ (*' + *0(MnY? m (k)j,(*A77 



We shall then define the primary source term operator as: 



S lm (k,v)= / d 2 ke~^ Vfa) 



e sw (Kv)-v b (k, v ) 



d_ 

drj 



(7) 



+ (V + &)(k, 7]) \ Y* TO (k) , (8) 



where we stress the fact that we have included the optical depth to Thomson scattering in the 
definition of the source term. In terms of Eq. ([8]), Eq. (|7|) can be written as: 



dV J (27r)V2 k S l™{k,r])jl(kAr]o) 



0) 
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In order to obtain a position-space description, we will express the coefficients Si m (k,r]) in terms 
of their counterparts in position space, by means of a Hankel transform: 

/dX 
-^^^X 2 S lm (X,v)jl(kX). (10) 

Using now the orthogonality relation for spherical Bessel functions, 

dkk 2 j L (ak)j L (bk) = - -j-^ 5 (a - b) , 

we obtain: 

e tl = r dr,S lm (X = A Vo , V ) . (11) 
Jo 

Eq. (jlip has a straightforward interpretation: it states that, in order for a source term at time 
r\ contribute to the CMB signal at time rj , that source must be located at the spherical shells of 
radius Arjo = rj — rj centered at the observer. The set of those spherical shells is a hypersurface 
which corresponds, of course, to the past light cone of the observer on {x ,?7 }. Since the visibility 
function is highly peaked at the time of decoupling, the primary source term contributes the most 
to the signal near the epochs when z{if) ~ 1100. 



B. CMB polarization in position space 



We shall now decompose the polarization as: 

Q + iU 



4/ 



(Vo, 6) = ^ Klrnijlo) 2 Y im (6) , (12) 
I in 

with the aim of determining the coefficients 7T; m (7/ ). The source terms in Eq. (jH are G 2 and a 2 . 
As a first step in an iterative solution, a 2 (which is of higher order in the visibility function) will not 
be taken into account - only the temperature quadrupole will contribute to generate polarization 
at this order. The first iteration is, therefore: 



X j2{kAm) JkAW > (13) 

where the time intervals are defined as Ar]\ = rji — r/ and ArjQ = r/ — r)\ . The source term Si m is 
the same as was defined in Eq. ([8j). Using once again the Hankel transform of Eq. (fT0|) we can 
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recast Eq. (|13p as: 

42 M 



'(1 + 2)! 



2^V (^- 2 ) ! 



'/<:. 



dm g(m 



j o Jo 



dfcfe ji(kX)—— — ^■j 2 {kAm 



(14) 



(fcA^) 2 

The interpretation of the expression above is the following: at time 77 a source term generates a 
temperature quadrupole. That quadrupole then generates, through a scattering at time 771 , the 
polarization which is finally observed at time r] Q . As we shall see, the integral of the second line of 
Eq. (|14p guarantees that the source term, at a distance X from the origin, is located in the past 
lightcone of the observer, for all possible m an d m The variable X will also be upper-bounded, and 
therefore the upper limit in the integration over the source terms will be replaced by a finite value 
that, as we shall see, corresponds the radius of the observer's past lightcone. 

The next step in the iterative solution is to take the n^l just computed and use it as a source 
term for the polarization itself - this means that now the polarization piece of the source term in 



Eq. ([4]), CK2, is no longer assumed to vanish. This contribution, which we will call 7rL:, 
given by: 

d 3 k 



is therefore 




dmoim) 



(2tt) 3 / 2 



dr)i g(m) i l — 



2 



a2{Km) 



Yf m (k) 



ji(kAr) ) 
(kAm) 2 



(2tt) 3 / 2 



where now the interval Am 
using Eq. (Ilip we obtain: 



rf 2 ke 2 (k,r ?2 )Y 1 ;(k) 



j 2 (k Am) 

(kA m ) 2 



24 

ji(k At?, 
(kAr] ) 



dm g(m 



m — V2- I n terms of the primary sources in position space, after 



'(1 + 2) 



>lo 



V a — 2)1 

dXX 2 S 2m (X,? ] ) 



dm g(m 



dm 9(in) / dr] 



9 / dkk 2 jiikX) 



ji(kAi] ) j2(kAr) 



1 J-32(kAr )2 ) 



(15) 



(kAm) 2 (kAm) 2 

At this point it is important to clarify our notation. We will always count the photon scatterings 
backwards in time: the time the photons are observed is always taken to be r/ D ; the last time that 
the photons scattered before being observed is 771; and so on. By convention, we will always evaluate 
the primary source term Si m at the instant 77, so the sequence of scatterings ends with 77. Hence if, 
as in the case described by Eq. (|15p . there are two scatterings between the generation of the signal 
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at r] and its observation at rj Q , then we have r] > rji > i]2 > r/, and the time intervals always express 
the differences between one time and the previous one, so in that case Atjq = T] — r]i, Arji = r]\—rj2, 
and A7/2 = f]2 — V- 

The same argument allows us to calculate the next order contribution, which takes into account 
three intermediate scatterings: 



r(3) ^ A 3 /(Z + 2)! 

'o Jo Jo 



rV3 r°o 

/ dq / dX X 2 S 2rn {X,rj) 
Jo Jo 



„ n2 f j, ; 2 • /. Y , 3i(kAvo) j2(kArn) j 2 (kArj 2 ) . , , . . 

x 9 J dkk 3i{kX) j2 {kA m) . (16) 

The general term in the iterative series expansion with n intermediate scatterings can be written 
as: 



(n) 3 /C + 2)! P° 1 f m 

^imiVo) = ~2^y (;_2)! y ^dim) x ^ _ j-yj y */2 ••• dr ?n T{ 5 ((7 ? 2) ■ - • £f(?7n)} 



/ dr? / dXJ 2 5 2m (X, 7? ) 
Jo Jo 



x / ^(«>f^f^ ... f^(^> (IT, 

* v ' 

(n— 1) times 

Here, T stands for the time-ordered product of the sub-intervals, whose purpose is to reproduce the 
chain of integrations mediated by visibility functions shown in, e.g., Eq. (|16D . 
The coefficients 7T; m appearing in Eq. (|12p can be expressed, therefore, as: 

oo 

vr im (r ?0 ) = ^vr i t ) (r ?0 ), (18) 

n=0 

with 7T^(rj ) given by Eq. flTTJ). 

An important condition for the validity of this perturbative expansion of the CMB temperature 
and polarization is that all the terms in the expansion of Eq. (|12p . with an arbitrary number n 
of intermediate scatterings, must be expressed in position space. However, this can only be true 
if the k integrals over products of spherical Bessel functions appearing in Eq. (|17p can be in fact 
performed, and are well-behaved. In the next Section we will show that, in fact, these integrals are 
probability densities for random flights with n steps, in a space of suitable dimensionality. 
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IV. RANDOM FLIGHTS AND THE CMB 



From the previous Section - specially from Eq. (fTTj) - it is evident that a complete treatment 
of the CMB in position space requires that some specific integrals of products of spherical Bessel 
functions should be computed. In order to fulfill this requirement, we will proceed in the following 
way: first, we will present a simplified version of those integrals, and we will show that they give 
rise to probability densities associated with random flights. Next, we will show how the integrals 
we have to solve can be expressed in terms of the random flight integrals. 

Let's recall the well-known identities satisfied by spherical Bessel functions: 



dz 



^J L+ |(z) 



Since the spherical Bessel functions are defined as: 



7T 




we are able to write: 



n-l 



dk k 2 j L (kr) J 



d 
dr 



L + 3 - 
2 



n—2 roo 



7T\(™+l)/2 



3L{kr n ) 



n— 1 r L+2 



kr\ L+ ^ 



~2~y 



[r(L + |)] 

n J 

(kr)l 



(kri) 



(19) 



The derivative of the second line in Eq. (fl~9|) is, in fact, the probability density associated with a 
random flight - see, e. g., Ref. Q. This is the probability density that a particle which moves with 
a constant (and finite) speed, and which starts from a given position in space, will be at a distance 
r from the point of origin, after changing randomly directions n times during its trajectory. The 
length of the intermediate steps are denoted by rj, i = 1, . . . , n. The order L of the spherical Bessel 
functions in these integrals is related to the dimensionality of the space where the flight takes place: 
namely, the dimension D of that space is given by D = 2L + 3. 
Following the notation employed by Ref. j(| we shall denote: 

n-l 



p n (r;ri, . . . , r n \ 2L + 3) := 



r [l + 




kr\ L+ ^ 



J L+ l{kri) 



(20) 



'o \ « J ~Ji (k ri ) LA 2 j 

However, the integral (1191) is still not what we need in order to solve the momentum integrals that 
appear in our iterations - see, e.g., the Fourier integration of Eq. (Qj 
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We will now extend the random flight integrals to include the scenario that appears in the 

context of the CMB. If I > L (which is always the case in our iterative solutions), then the product 

of two spherical Bessel functions of order I can be written in terms of a single spherical Bessel 

function of order L. Th.s is a conscience of GegentWs .elation fl Q and of the orthogonality 

of associated Legendre polynomials - see [5| for a derivation: 

(_-\\L rX+R / X E\ L ~ 1 

ji^UkR) = ±-L- / drk L Pf L (cosa) (sina) L j L {kr) (21) 

2 J\x-R\ \ r J 

where r, R and X must form a triangle, with the angle a being given implicitly in terms of the 

relation r 2 = R 2 + X 2 -2RX cos a, and (cos a) is an associated Legendre polynomial. Applying 

Eq. (|21|) and Eq. (|19p with L = 2 we can recast the Fourier integral in Eq. (|17|) as: 

F l2 ^ m ,x-^...Ann) ■= ^^^^^^f^nfi^^^) 

= r / ci(cosa) I 5-^ ) P,~ 2 (cosa) sin 2 a 

x p n (r;A? ? i,...,Ar ?n |7) , (22) 



where we used the notation introduced in Eq. (I20p . 

In Eq. (|22p we have introduced the function F^iArjo, X; Ar]i, . . . , Aij n ), which denotes what we 
shall call extended random flight integrals. In general, Fil is defined by Eq. (lAip . 

We should examine Eq. (I22p more carefully. As anticipated, the presence of 
p n (r; Arji, . . . , Ai] n \ 7) should not be surprising, due to the interpretation of a random flight 
process and its validity with respect to the physics of the recombination. The dimensionality 
(D = 2 x 2 + 3 = 7) of the space associated with the random flight, however, is not yet fully 
understood. That dimension is determined by the order of the spherical Bessel functions which 
mediate the sources of anisotropics and the final CMB signal, but since only the quadrupole of the 
temperature fluctuation contributes to the polarization, the spherical Bessel function of order 2 is 
the one that characterizes the random flight for the CMB. A possible explanation for this dimension 
is that, after separating the angular dependence of the CMB from its radial and time dependence 
through the spherical harmonic decomposition, the light cone has only two dimensions left. Since 
the multipole L has 2L + 1 degrees of freedom, we end up with 2L + 3 dimensions where our relevant 
variables can perform random flights. However, a more refined argument to explain the dimension 
7 is not yet known. 

Looking back now at Eq. (I17p . we recognize that for the 7T,", term in the polarization expansion, 
the probability density p n (r; Ar)x, . . . , Aij n \ 7) appears clearly. The interpretation of an expansion 
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in the number of interactions during recombination is therefore strengthened. We should remark, 
however, that explicit formulae for p n (r; Arji, . . . , Arj n | 7) are not known in the general case 

As discussed above, the intervals A771, . . . , Ar/ n express the times elapsed between consecutive 
scatterings. All these subintervals are elements of a partition of the time interval rj\ — r], which is 
therefore the total time elapsed since the instant the photon leaves equilibrium with matter, at time 
rj, until the instant r/i when the photon has last scattered prior to its observation. This time interval 
represents, therefore, the effective duration of recombination for a given photon. The lengths of 
each subinterval are weighted by the visibility functions, and integrated in order to contemplate all 
possible histories for photons during recombination. The interval r/ — r/i expresses the time elapsed 
since the photon scattered for the last time, before it is observed at the time r\ (we are assuming 
that no further scatterings take place during this interval). We can represent a photon's history by 
means of Fig. [TJ on the left we show the photon's interactions prior to observation at the vertex 
of the cone, while on the right we show a diagrammatic representation of that history, with the 
relevant elements that appear in Eq. ()22D . 

We call special attention to the diagram represented on the right of Fig. [TJ This diagram 
represents the extended random flight performed by a photon during recombination. The steps Arj n , 
. . . , A771 (going forward in time) belong to a standard random flight, and describe the trajectory 
of a photon that has left equilibrium with matter at an instant i], then propagated freely for 
a distance (or a time interval) Ar/ ni then Thomson-scattered with an electron at time r/ n , then 
propagated freely for a distance Ar] n ^i, and so on until the instant 771, when it scattered for the 
last time. The standard random flight ends at the instant r/i. The photon, at that moment, is 
a radius r away from the point where the flight started. The steps indicated by A770 and X do 
not belong to the standard random flight, but are present in the extended random flight, and are 
introduced through the spherical Bessel functions of different order in the k integral of Eq. (j22|) . 
This is necessary because those two steps are not associated with any movement between successive 
scatterings, but are in fact associated with the distance from the observer to the origin of the 
photon, and to the end-point of the random flight. It should be strengthened that Ar/o, X and r 
are related by r 2 = At]q + X 2 — 2 At]q X cos a. Since < r < Ar/i + A 7/2 + . . . + Ar/ n , it follows 
that < X < Ar/o + Ar/i + Ai]2 + . . . + Ar/ n , which then determines the domain of dependence of 
the problem. 
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Figure 1: Connection between scatterings of a photon during recombination (left) with a diagrammatic represen- 
tation of the random flight (right). The steps Arj„, . . . , Arji correspond to the lengths of the photon's trajectory 
between successive scatterings. A770 corresponds to the propagation since the photon's last scattering during recom- 
bination, at time 171 , and the observation at time r\ . In both diagrams the observation takes place at the upper 
vertex. 

V. THE POLARIZATION IN POSITION SPACE 

We shall now go back to Eq. (|17|) . In terms of the extended random flight just introduced, the 
polarization coefficient to n th order, ^^(Vo), can be written as: 



<„) 3 (l + 2)! [«" 1 rm 



n irn = -J- \ {] _ 9V / dr li 9(m) — TTT / dr l2 ■ ■ ■ drj n T{g(r]2) . . . g(vn)} 



nn-l /7r N \(«+l)/2 rrjo-rj 

[r (I)] Jo Jo 

d(cosa) ( 3~~ J ^z _2 ( cos a ) s ™ 2 a Pn( r i ^7/1, ■ ■ ■ j Ar?n I 7) , 

where we have already used the aforementioned upper bound for the variable X. 
Eq. (|23p can be understood as the combination of three procedures: 



(23) 
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• The integration over a corresponds to a marginalization over all possible paths composed of 
n steps of lengths Ar/i + . . . + Ar/ n that have a net displacement r determined by X and A 770, 
as shown in Fig. [2j This "average over paths" is a function of X, A770, Ar/i, . . . , Ar? n . 




Figure 2: Marginalization over all paths with n steps, composed of the intermediate displacements A771, . . . , Ar/ n , 
which lead to a fixed displacement r with respect to the origin of the flight. The distance r is determined by X and 
A770 for all possible angles a. 

• The contribution from the quadrupole of the primary source term, S2 m (X, 77), is then medi- 
ated by this "average over paths" that was just described, for all possible values of X. The 
maximum value that X may reach, is Atjq -\-At}i ~\~ . . . H~ A?y n = 7? — 77 , which is nothing but the 
radius of the observer's past light cone up to the time 77. After computing the contribution 
of the source terms, we end up with an expression that is a function of A770, A771, . . . , A?7 n . 

• The last step is to let the intervals A770 , A771 , . . . , At7„ assume any values through the inte- 
grations, each one weighted by its corresponding factor of the visibility function to take into 
account the probability that the photon will scatter at that instant of time. This accom- 
plishes the goal of accounting for the contribution from sources at all distances, and over any 
possible number of intermediate steps of the extended random flights. 
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Adding the contributions from all 7r," we obtain: 



Kir 



3 (i + 2)\ n° f l r 

gn— 1 /7r\ (n+l)/2 „,j n 



d7?2 • ■ • rfr/ n T{g(r] 2 ) . . . g(Vn)} 



[r(D] 



7M n— 1 



./Q 

2 



d(cos a) ( g-^ j P^ _2 (cos a) sin 2 a p n (r; A771, . . . , Ai] n | 7) . 



In terms of the extended random flight integrals, Eq. (J22J), we can write: 



3 /(Z + 2)! 



2tt y (Z-2)! J '^'"^ (n-1)! Jo 

/•J7„ roc 

/ d7/ / dXX 2 5 2m (X,? ? )F i2 (Ar ? o,X;Ar ? i,...,A^ 
Jo Jo 



(24) 



00 on— l rvi 

dm g(m) E 7— — I d m ... d^iVfaV^ . . . 5 (r? n )} 



(25) 



VI. DIAGRAMS 



We are now able to introduce a diagrammatic representation for all the terms that appear in 
the series expansion of the temperature and polarization at n-th order. The graphic elements are: 

• Solid lines mean polarization, and dashed lines mean temperature; 

• The vertical lines to the right of the diagrams represent the observable being calculated (time 
runs upward); 

• The interception of two lines determine what sources are being considered for a given observ- 
able. 

The rules for constructing these diagrams will be detailed in the next Section, but it is already easy 
to represent diagrammatically the contribution to the CMB polarization computed in Eq. {\24\i : 



-9 + 



+ 



-9 + 



+ 



(26) 
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VII. FULL COUPLED TEMPERATURE EVOLUTION 



We will now show that the formalism developed in the previous Sectons is not spoiled when all 
the contributions to the temperature fluctuations are included. Let's go back to Eqs. ([I]) and ([3]), 
and write: 



>lm 



2i l dri , 

/o J (2vr)V 



dk 



d 2 



e sw (k, V )-V b (k,r,)^- 



+ e -^)(*' + $')(k,r ? ) 



1 + 3- 



d{kAr, ) 2 _ 

The first two terms in this equation have already been treated in Section III. We now pro- 
ceed to analyzing the remaining term, i. e., that which contains the source P(k, rf) = 
\ [02 (k, rf) — \/6 «2(k, ?7)] . We write, therefore: 



'Ira 



JO 

oo 

□ (0) , \^ n(n) 
7 lm "i" lm ' 

n=l 



172 ^ / d'kg( V )P(k,r,)Yt m (k) 



1 + 3 



2 



d{kAr, ) 2 



ji(kArj ) 



(27) 



where we have defined the 9) m in terms of the integral over sources. 



A. Temperature at first order 



We shall now calculate the lowest-order contribution described in Eq. (|27|) by taking P(k, rf) — > 
^Q 2 (k, if) - i. e., only the temperature quadrupole is taken as a source for the temperature fluctua- 
tions, since the polarization (which is of higher order, since its source is the temperature quadrupole) 
is at least a second-order contribution to the temperature. The first order contribution to the tem- 
perature, is therefore computed as: 



7 lm 
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u dr?i5(7?i) y (^ 
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ji(kAr) ) 
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dm g(m 



rvi r°o 

/ d V / dXX 2 S 2m (X, V ) 
Jo Jo 
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1 + 3: 



t) 2 



ji(kAri ) j 2 (kAr]i 



(28) 



'd(kA Vo ) 2 

where we have once again transformed coefficients into the position-space description using Hankel 
transformation. We must now treat the integral in the last line of the equation above: 

d 2 



H2 



dk k 2 j t (kX) 



1 + 3 



d(kAr] 



ji(kAri ) j 2 (kArn) . 
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Clearly, we can write: 



112 



o 2 







9(Ar ?0 )2 Ar] d(A Vo ) 



(A Vo y / dkk'j t (kX) 



ji(kAr) ) 
{kA Vo y 



hikAr/i 



+ (A Vo y / dkk^MkX) 



ji{kArj ) 
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j 2 (kAm 



(29) 



which we can then express as: 

d 2 1 d 



In 



d(Ai] ) 2 Ar] dAr] 



[(Ar ?0 ) 2 F, 2 (Ar ?0 ,X;Ar ?1 )] + (At? ) 2 G$ {A Vo , X; A m ) , (30) 



(2) 

where the integral corresponding to G l2 ls related to the already introduced function F\2 - see Eq. 
(|22p . In fact, using an identity that is proved in the Appendix, Eq. (|A7|) . we obtain: 



9 (l) 

7 lm 



i r h 

2^ 



I drj ig ( m ) dr, dX X 2 S 2m (X, V ) 
o Jo Jo 

[(Ar, ) 2 F l2 (Ar, 0) X;Am)] 



2 



1 d 



d(A Vo ) 2 Ar, dAr, 



+ 



1 











+ 



(A Vo X) l+2 F {l+1)2 (A Vo ,X;A m ) 



(31) 



(Ar] ) l X l + 2 dX \dAi] Ar, 0/ 

In conclusion, first-order corrections to the temperature can also be written in terms the probability 
densities of extended random flights. 



B. Temperature at second order 



If we want to calculate further contributions to the temperature, then we must include the 
contributions to that come from the temperature quadrupole, given in Eq. (|30p . and those from 
«2> given in Eq. (fH|) . Combining those two contributions leads to: 



ft) 

'Ira 



dm aim) / dmgim) / dr] dx x 2 s 2m (x,r]) 

Jo Jo Jo 

d 2 

1 + 3- 



1 

^ JO 

x j dkk 2 j l (kX)j 2 (kAr ]2 

+ -/ dmg(m) d mg (r, 2 ) dr, dx x 2 s 2m (x, v ) 

7T Jo Jo Jo Jo 



j 2 {kAr] 



x / dkk 2 j l (kX)j 2 (kAr ]2 [ 



d(kAm) 2 
) 

(kAm) 2 f ' ~d(kAr, ) 



1 + 3 



d 2 



d(kAr] ) 2 



ji(kAi] 



1 + 3 



<9 2 



ji(kAr] ) . 



(32) 



(2) 

The first term comes from the contribution of 02 to 9)' , and the second term takes into account the 

(2) 

contribution of a 2 to 9)' . These contributions can also be expressed in terms of the extended ran- 
dom flight integrals: terms with two derivatives can be expressed in terms of Fi 2 (Arjo, X; Am, Ar/ 2 ); 
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terms with one derivative can be expressed in terms of 2(A?7o! X; Ar/i, A772); and terms with 

no derivatives can be expressed in terms of -F^+2) 2(Ar/o, X; Ar/i, A772), as indicated in Sec. IA 21 

The complete series expansions for the temperature coefficients, 9i m , is represented diagram- 
matically, up to second order corrections, as: 



9(2) 



+ 



I 
I 
I 
I 
I 

-6 

1 
1 
1 
1 
1 



+ 



1 
1 

-6 
1 
1 
1 

^> 
1 
1 
1 
1 



+ 



1 
1 
1 
1 
1 

-6 
1 
1 
1 
1 
1 



(33) 



The first diagram represents 0$, the second, 9\t^ - see Eq. f[30|) - and the two last diagrams 
represent 9}', as shown in Eq. (|32p . 

Even if we do not calculate here explicitly the contributions at all orders, the results presented 
in Sec. IA 21 show that all corrections can be expressed in terms of extended random flight integrals. 

-(n) 

In general, 9) _ will be given by: 



(i) 



n(n) 
7 lm 



J2 0fi( Ar 7°' • • • ' Ar >n,X) F {l+qp (Ar] , X;A Vl , 



A?7 n ) , 



q=0 



where we represent by 9^ the coefficients that appear in the expansion. As was seen explicitly in 
the case n = 1 , these coefficients are differential operators combined with powers of the subintervals 
of (ijo^v) an d of X. The complete expression for 9^ will also depend on knowledge of the expression 
for polarization up to the order n — 1 - which is the reason for leaving the superscript n explicit in 
9\ n) . 

Lmq 



C. Polarization at second order 



The last step in our construction of the series expansions of the CMB in terms of extended 
random flight probabilities is to treat the contribution of the temperature quadrupole to the polar- 
ization, allowing for the corrections to the temperature that come from higher order contributions 
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of the source terms. We shall write, then: 



7T 
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k 2 S 2m {k,r]) j 2 (kAr]i) 



(2tt)V2 
j 2 {kAr]i) 



k z / ^kP(k,r? 2 )YL(k) 



1 Jf (fcAr/p 



(34) 



where the bar over 7rL: means that we are taking all corrections up to second order. 

At this point we can separate the first-order and the second-order terms, and focus on the 
contribution from Q 2 to P: 
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j 2 (kAr)i) 
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1 + 3 



d{kArji 
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(kAriof 
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jijkArjo) 
(kAri ) 2 



1 + 3 
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d(kAr/i 



j 2 (kAr]x) j 2 (kArj2) > . 



(35) 



The last integral of the equation above, between curly brackets, can be expressed, after using Eq. 
(fA7j) . as follows: 



j dkk 2 n (kX) 



ji(kAr] ) 
ikA m ) 2 
(Arn) 2 d ( d 



1 + 3 



+ 



d 2 



d(kAr]i 



j 2 (kAr]i) j 2 {kAr} 2 ) 



iAr] X) l + 2 dX \dArjo Ar] 

, * i a 



(Ar? Xy+ 2 F {l+1)2 (Ario,X; A m , A m 



iA m y F l2 (Ari ,X;A m ,Ari 2 ) 



(36) 



\diAr]i) 2 Ar]idAr]\ 

As we can see, despite the somewhat cumbersome coefficients, the final result can be expressed, 
again, as a combination of the extended random flight functions Fi 2 (At]q, X; Ar/i, Arj 2 ) and 
F (i+i)2i A Vo, X; Arix, Arj 2 ). 
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The last term in Eq. (I35p has the same order as it)' . Therefore, to second order in temperature 
corrections, the polarization results from the following three diagrams: 



(2) 
Im 



+ 



+ 



—9—6 



(37) 



As we showed above, this contribution can be expressed fully in terms of extended random flight 
probabilities, therefore, up to second order the polarization coefficients can be written as: 



*Zm = 7r /mo(Ar/o, Ar/i, Ar] 2 , X) F l2 (Ar] , X; Ar/i, Arj 2 ) 

+ n lml (Ari , Arn, A?] 2 ; X) F( l+1)2 (Ai] , X; Arn, Ar] 2 ) ■ 



(38) 



The coefficients Tri mq (q = 0, 1) are differential operators, and depend on the combinations of the 
subintervals to some power. The coefficient ni m o takes into account the contributions from the 7T^ 
term appearing in Eq. (|35p and also the contribution comming from the second term in the right 
hand side of Eq. (|36|) . The coefficient ixi m i can be obtained by working out the first term in the 
right hand side of Eq. (|36|) . 



If we examine Eqs. (f35|) and (|37p . it can be seen that 7fP^ is given by the sum of and the 
contribution from the temperature corrected by the polarization. We can write Eq. (I37p in a more 
symmetric way if we note that the following diagrams are equivalent: 



(39) 



This equivalence is due to the fact that both diagrams have the same vertical lines at the right (we 
are calculating contribution to the polarization); they have the same number of vertices; they have 
the same kinds of lines connecting the vertices (continuous lines); and they have the same initial 
source, i. e., the temperature fluctuations (the dashed lines). Hence, it is better to write Eq. (|37l) 
in terms of the diagram on the left-hand side of Eq. (139 j) . since it is easier to draw these types of 
diagrams when there is a large number of scatterings, as shown in Eq. (f26|) . 
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In general, to order n: 
n-1 

= E 7r £( A %, Ar?!, . . . , Ar? n , X) F (l+q)2 (A m , X; A m , . . . , A ??n ) . (40) 

9=0 

(n) 

Explicit formulae for all the coefficients 7Tj m are not known explicitly, but can be computed iter- 
atively after the identification of the diagrams that contribute to some order, and by employing 
the techniques presented in Sec. IA 21 Naturally, the computation of 7f&' requires the knowledge of 
9^ , the computation of ^ requires the knowledge of both 2 "* and 7f[™ , and so on and 
so forth - therefore, in this sense these coefficients form a hierarchy. 

D. Diagrammatic rules 

If we accept that Boltzmann's equations, as presented in Eqs. (JTJ) , ([3]) , d2J) an d ©, are accurate 
enough to allow the computation of the CMB temperature and polarization to all orders, then we 
can ask how the diagrams for our series expansion should be built. Clearly, when further precision 
is required, quadratic and even higher-order corrections for the metric perturbations may have to 
be included in Boltzmann's equation, and at some point the diagrams presented here may no longer 
improve the accuracy of the result. However, if we assume that the Einstein-Boltzmann system of 
equations is exact, then all contributions can be written as a sequence diagrams, whose building 
blocks are presented below: 
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diag. segment numerical factor 



Bessel function 



_3_ / {1+2)1 
'2tt V (1-2)1 
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1 4- q 9 2 .. 
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n n q J2(fcA) 

i 
i 
i 

As a general rule of thumb, the functions ji(kX) and j2(kAr/ n ) must appear as the first and 
last elements inside the k integral. Here, r] n stands for the last index in the partition of the time 
interval (r/ D — r/). The intervals A must be replaced by the convenient difference of consecutive 
conformal times which represent the initial and the final instants of that segment of the random 
flight. Outside the k integral there should be the integral over the spatial position of the primary 
source, J dX X 2 S , 2 m (X, rj). Finally, we have the integrals over conformal time weighted by visibility 
functions - except the integral over T], which does not carry a factor of the visibility function. 

For example, the following diagram represents the third-order contribution to polarization that 
arises from the polarization at second-order, when that second-order polarization is due to the 



22 



first-order correction to the temperature: 
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(41) 



The diagram below, on the other hand, represents a third-order contribution to polarization 
arising from the second-order correction to the temperature, when that second-order correction 
itself comes from the first-order correction to the temperature: 
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i i 
i i 
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x / dXX z S 2m (X, v ) { dkk'MkX) 



ji(k Ar] 
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VIII. FOURIER-BESSEL SERIES FOR THE EXTENDED RANDOM FLIGHT 

PROBABILITIES 

From the previous arguments it is obvious that a full description of CMB temperature and 
polarization depends upon the evaluation of the extended random flight integrals F[ 2 . As we 
showed in Eq. (|22~1) . the extended random flights can be expressed as marginalizations over random 
flight probability densities. Unfortunately, explicit formulae for random flight probability densities 
in seven dimensions are not known, which means that at orders higher than three the computations 
start to become highly complex. As a computational tool to deal with the extended random flight 
integral, we will introduce their series expansion in terms of a Fourier-Bessel series - this was also 
suggested by [jjj]. 

We will expand the function Fn,(R, X; r±, . . . , r n ), which is ubiquitous in our previous discussions 
(just substitute Aryo —> R, Ar/i — > r±, etc.), and which was defined in Eq. (JATJ) , in terms of a Fourier- 
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Bessel series. Given the behavior of this function and its interpretation, the segments R, n, . . . , r r , 



are fixed for any particular flight. Since r = \/R 2 + X 2 — 2 R X cos a < t\ + + . • • + r n =: S n , it 
follows that the variable X must have itself an upper limit - see Fig. 1 for the geometry involved in 
these limits. In fact, in the limit where all flight steps r\, . . . ,r n are collinear, we have that r = S n , 
and in that case X, R and S n form the sides of a triangle. Hence, the maximum value of X is 
determined jointly by both R and S n , as < X < R + r\ + T2 + • • • + r n := S n . 

Now, we know that well-behaved functions with a limited domain can be expressed in terms of 
a Fourier-Bessel series fisl . In fact, since Fn, is one such well-behaved function, we have: 



F lL (R,X;r u ...,r n ) = J2tiujl(4 ] x) , (43) 



y( n ) „\ ( J n ) 

i=i 

» - J id a „A J 



where x Li = zySn, and z\ is the i-th root of the spherical Bessel function ji(z). The coefficients 

flu are s iven b y 



fiLi = 03 -2 f i, ["dXX 2 3i (4° X) F lL (R, X- n , . . . , r n ) 
Using Eq. (f22|) we obtain: 



(44) 



i=l b nJl + M) f x [f R ) q=l f x W Tq 



(45) 



Now, Eq. (I2ip allows us to write (j45J) in the form: 

(—1) L 7r f' Sn /RX\ L l 1 
F lL (R,X; ri ,...,r n ) = J q dr P" L (cos a) (sin a) L 



x 



1=1 J l+1\ I' g=l lj ( v . ' r q \ 

where r 2 = i? 2 + X 2 - 2 ii X cos a. 

The extended random flight integrals can, therefore, be much more easily computed with the 
help of Eq. (|46p . since in that representation only a discrete number of modes need to be added, 
instead of the full-fledged integral that defines the random flight probability density. 



IX. LARGE NUMBER OF SCATTERINGS AND THE H-THEOREM 

In this Section we shall address a more fundamental question. Up to now we have presented a 
formalism that allows us to express the CMB polarization and temperature in position space, in 
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terms of an expansion over the number of interactions the photons have suffered during recombina- 
tion. However, at this point we can ask: up to what order the contributions should to be taken into 
account in order to yield an accurate expression for the physical observables? But this question 
can only be answered if a concrete problem is given, so there is no general answer. Hence, we will 
invert this question, and ask: how important is it, for the observables, that we reach a certain order 
n in the expansion? To be more specific, we will study the contributions arising from very high n 
terms in Eq. (|24D . 

Since all computations depend on the probability densities for random flights with n steps, we 
will analyze the behavior of p n { r ] A7/1, . . . , Ar/ n | 7), for large n. Although explicit formulae for these 
probability functions are not known, an asymptotic expression is known in the case where all steps 
have the same lengths, Ar]i = Ar/2 = . . . = Ar/ n =: A, and n — > 00 Q|. In that case we have: 



is the confluent hypergeometric function [151 ] . Expanding \F± in a power series we obtain: 




(47) 



where: 





(48) 



Now, consider the sum: 




It follows, therefore, that: 




(49) 



Consequently, for fixed r > 0, A > 0, and r < nA 



lim p n (r;A,...,A|7) =0. 
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Since p n is a normalized probability density, we conclude that, in the limit above, the probability 
distribution collapses into a Dirac delta-function centered at the origin, which is in agreement with 
the general features attributed to the random flight probability densities In other words, in- 
creasing the number of changes of directions (i.e., scatterings) in a random flight leads to trajectories 
with increasingly smaller displacements from the origin. 

Let's investigate the consequences of this fact for the polarization through Eq. (|23p . Consider, 
in that respect, the integral: 

2 
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n— >oo 
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(50) 



We know, however, from either l6fl or l5fl, that: 



and also that 



dk ji(kri) ji(kr 2 ) j2(kr 3 ) = - ^r- P, (cosa)(sina) , (51) 

4 Tg 



roo r 2 r 2 ^ 

lim / dk jiikri) ji(kr 2 ) j 2 (kr 3 ) = -f -f -S(n - r 2 ) . (52) 
J r\ 15 2 

Therefore, we have that: 

5 2m (r?o,r?) = £° V dXX 2 S 2m (X,r 1 )^^5(Ar ]0 -X) 

= voj^r S2m{Ar]o ^- (53) 

The final lesson is that, in the limit of an infinite number of scatterings, we basically end 
with instantaneous recombination - with the crucial caveat that the source term is damped 
by (A/A?7o) 2 — > 0. But in instantaneous recombination the temperature fluctuations have no 
quadrupole term which means that the source for the polarization vanishes. Similarly, for the 

(n) 

temperature fluctuations the derivatives of the integral which appears in 9)' ensure that the signal 
vanishes as well. 

The large n limit, therefore, cannot contribute significantly to the final observables. If we regard 
this order as the number of scatterings that a photon has experienced during recombination, then 
the collapse of the probability density associated with random flights expresses a fundamental 
physical fact: increasing the number of scatterings takes the system closer to thermodynamical 
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equilibrium. As we know, in thermodynamical equilibrium the CMB temperature is given simply 
by the Planck distribution, with no distortions (apart from the dipole which is induced by the 
velocity of baryons). Maximizing the von Neumann entropy also ensures that, in equilibrium, no 
net polarization is generated. These two facts, together, constitute the expression of Boltzmann's 
H-theorem to the CMB Q- 

Describing the CMB by means of random flight probability density functions provides, therefore, 
an illustration of the statement of the H-theorem: when few scatterings take place during recombi- 
nation (low n), photons and electrons are surely out of equilibrium, and each interaction generates 
temperature fluctuations and polarization that are not sufficiently erased by the subsequent scat- 
terings. As a consequence, the low-n terms in the expansions account for the largest contributions 
for the signal of the physical observables. For a large number os scatterings, the signal is washed 
out by further scatterings, which can be expressed in terms of the collapse of the probability density 
in the large- n limit. 

All that is left to discuss is the physical grounds for assuming (as we did in this Section) that all 
the steps performed in the random flights have the same length. Despite the fact that this working 
hypothesis was not proposed because of any physical reason, but purely because of mathematical 
convenience, we can nevertheless argue that, when the electron- photon system was close to the 
tight-coupling regime, the assumption that the visibility function is step-shaped is not completely 
crude. Continuity with respect to the arguments of the probability density distribution makes our 
working hypothesis less unnatural. 



X. CONCLUSIONS 



The Boltzmann hierarchy for the problem of the evolution of temperature and polarization 
fluctuations of the CMB is equivalent to a system of coupled integral equations. We showed that 
this system can be written in position space, and that the objects that appear in this description are 
the same probability density functions that appear in random flight problems. The emergence of 
random flight probability functions from Boltzmann's equations clarifies the physical interpretation 
that CMB temperature fluctuations and polarization are generated from scatterings of photons by 
low-energy electrons during the recombination. More importantly, we showed that the number of 
scatterings during recombination is a key ingredient in the description of this problem - in fact, a 
perturbative expansion can be performed in terms of the number of scatterings. 

We showed, using asymptotic formulae for the random flight probability distribution, that con- 
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tributions coming from high-order terms (i.e., many scatterings) should indeed be negligible. Since 
high-order terms represent the past history of photons that have scattered many times during re- 
combination, we concluded that their vanishing contribution to the CMB signal is an illustration 
of Boltzmann's H-theorem. In the context of the CMB, this theorem states that, in thermodynam- 
ical equilibrium, the CMB temperature should be given by a Planck spectrum, and that the net 
polarization should vanish. 

Using the fact that the random flight distribution functions vanish identically if one cannot form 
a closed polygon from its intermediate steps, we presented Fourier-Bessel series expansions for the 
associated probability distribution function. These expressions lead to simple numerical recipes for 
the computation of these distributions, since explicit formulae for them are not presently known for 
any dimension. Another very important result which derives from the emergence of random flights 
is the fact that, at each given time, the domains of dependence of the problem are compact sets. 
This is not obvious from the usual treatment of Boltzmann's equations in Fourier space. 

From a more formal perspective, it would be important to understand which classes of Boltz- 
mann's equations are amenable to a treatment in terms of random flight probability densities. Many 
of those problems could be then examined in light of the formalism that has been developed for the 
study of the CMB. 

Finally, we can foresee some possible applications of this work. The series expansion in terms of 
the number of scatterings can be used for numerical simulations of constrained maps of temperature 
and polarization. It may be sufficient to consider only up to three or four scatterings of the photons 
in order to achieve an accuracy which is sufficient for most applications. The partial summation of 
some classes of diagrams (in the spirit of a renormalization group) can also be performed, which 
would lead to more accurate results without increasing the complexity. In particular, the position- 
space analysis is especially suited for introducing gravitational lensing (work in progress). 



We shall show the possibility of expressing a larger set of integrals in term of random flight 
probability densities. Specifically, we will concentrate in expressions with the form: 



Appendix A: Extended random flight integrals times even functions 



/ 

Jo 



oo 



ji(kR) 
(kR) L 



n-l 



3L{kr q ) . 
(kr q ) L 3 ' 



F lL (R,X; ri ,...,r, 



) 



dk k 2 ji{kX) 



n 



(k 




(Al) 



? =1 



28 



Our aim is to express integrals like: 

G lL (R, X; n, . . . , r n ) = J™ dk k 2 f(k) ji(kX) ^| fl J^-j^krn) (A2) 

in terms of combinations of Fy^^R, X; n, . . . , r n ), for some adequate I' and L'. In Gil, f(k) is 
taken to be a real function such that its Taylor series carries only terms with even powers of k. The 
simplest such function is ak°, for a G R. Naturally, in that case the connection we are looking for 
is quite obvious. We will, therefore, study monomials in the Taylor expansion, namely, we will take 
f(k)=k 2 . 



1. The case of Gil{R, X; r\, . . . , r n ) when f(k) = k 2 



The following will be the object of our attention: 
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where the superscript indicates that we are treating the power k in the expansion of some f(k). 
Recall the recurrence relations obeyed by spherical Bessel functions: 

^-[x» +1 j u (x)] =x» +1 j^ 1 (x). (A4) 

It then follows that, for some integer that: 



x l+2 Mx) = ±( x l + 2 j l+1 (x)) . 
ax 



The chain rule gives: 
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Inserting ()A5P into ()A3P we obtain: 
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However, 
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for any function u(R). Hence, 
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We obtained, therefore, a relation: 
Gg } (Ji,X;ri,...,r n ) 
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(RX) l+2 F {l+1)L (R,X; ri ,...,r n ) , (A7) 



(RX) l + 2 dX \dR R / 

i. e., G\^'(R, X; r\, . . . , r n ) can be expressed in terms of derivatives of Fn + ^ L (R, X; r\, . . . , r n ), 
which is the result we were looking for. 



2. Higher powers of f(k) 



We will not derive a closed-form expression for G[{(R,X ;ri,...,r n ) in the case of any even q, 
however, we will illustrate how to generalize the argument presented in the previous subsection to 
these higher powers. 

First, we point out that higher powers in the expansion of f(k) must be canceled out by corre- 
sponding factors 1/k when Eq. (jA5]) is iterated for the appropriate number of times. For the fourth 
degree monomial in the expansion of f(k) we would have: 
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This holds for ji(kR), as well as for ji(kX). We must bear in mind that, in the original equations, 

terms like appear, and therefore we must change the orders of the derivatives with respect 

to R with a term like (kR)~ L . This can be performed easily with the aid of the Leibniz formula 

for the re-th derivative of products of functions: 
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This process must be iterated until all derivatives with respect to R have been interchanged with 
terms like (kR)~ L . 

As a final remark, we note that, even if the analytical problem is quite cumber- 
some, G^f L (R, X;ri, . . . ,r n ) can always be obtained, for even q, in terms of derivatives of 
F (i+ q /2)L(R,X;ri,. . . ,r n ). 
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